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TIME REVERSAL FOR RADIATIVE TRANSPORT WITH 
APPLICATIONS TO INVERSE AND CONTROL PROBLEMS 

SEBASTIAN ACOSTA 

Abstract. In this paper we develop a time reversal method for the radiative transport 
equation to solve two problems: an inverse problem for the recovery of an initial condition 
f^ ' from boundary measurements, and the exact boundary controllability of the transport field 

CN . with finite steering time. Absorbing and scattering efi'ects, modeled by coefficients with 

$H ' low regularity, are incorporated in the formulation of these problems. This time reversal 

Mh, approach leads to a convergent iterative procedure to reconstruct the initial condition pro- 

■^ ' vided that the scattering coefficient is sufficiently small in the L°° norm. Then, using duality 

("^ . arguments, we show that the solvability of the inverse problem leads to exact controllability 

of the transport field with minimum-norm control obtained constructively. The solution 
approach to both of these problems may have medical applications in areas such as optical 
^ ' imaging and optimization of radiation delivery for cancer therapy. 

< 

c^ ; 1. Introduction 

The radiative transport equation, also known as the hnear Boltzmann equation, gov- 
erns the propagation of particles as they interact with the underlying medium. Hence, this 
^ \ equation has applications in various scientific disciplines, including gas dynamics, optics, 

r~;~ ■ astrophysics, nanotechnology, and biology. The mathematical modeling for some of these 

applications is developed in the following books [H El El IH [S] among others. This publica- 
cn i tion is primarily motivated by the connection between transport phenomena, and medical 

t:;}^ \ applications such as optical imaging and the optimization of radiation therapy. More pre- 

O ' cisely, we construct a time reversal method to solve two problems governed by the radiative 

transport equation. These problems are the recovery of an initial condition from boundary 
measurements, and the exact boundary controllability of the transport field. We shall refer 
to the former as the inverse source problem, and to the latter as the control problem. 



cn 



'i 



1.1. The inverse source problem. The understanding of transport problems leads to 
imaging methods for biological media with particular applications in medicine [SI El |B1 13 HDl 
[TT] . Another area of interest is the development of wave-based imaging techniques suited 
for turbid or heavily cluttered media modeled with stochastic differential equations. The 
connection to the radiative transport equation is provided by the fact that the Wigner 
transform of stochastic wave fields satisfies the transport equation in specific asymptotic 
regimes. For details and recent developments in this theory, we refer to ^2\ [T3l [TH [151 UHl 

[ni[18l[l9lE0l[2llE2]. 

The general radiative transport problem is governed by the following system, 

if + (^-V)M + /iaM + /is(/-/C)M = / in (0,oo) X (fixS), 
u = uq on {t = 0} X (i7 X §), 

augmented by a prescribed in-fiow profile which we assume to be vanishing. The properties of 
the medium are modeled by the absorption coefficient /la, the scattering coefficient /Xg and the 
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scattering operator /C. The precise definitions and assumptions concerning these coefficients 
are made in Section [2l The transport equation has solutions u = u(t,x,9) representing the 
density of radiation at time t G [0, oo), position x E fi G M", moving in the direction 6* G S at 
speed c > 0. Here S denotes the unit sphere in R"^. The driving sources in this problem are 
the actual forcing term in the right-hand side of the equation denoted by / and the initial 
condition uq. 

For inverse source problems in transport theory, the goal is to reconstruct some of the 
driving sources in the formulation of the transport problem from the knowledge of out- 
flowing boundary data. This theory can be divided into two main areas, namely, transient 
and stationary problems. To the best of our knowledge, there is no rigorous work concern- 
ing the reconstruction of the initial condition uq from boundary measurements for general 
heterogeneous, absorbing, scattering media. In this paper, we partially ffil that void by 
developing an iterative method, suited for weakly scattering media, to recover the unknown 
initial condition. 

Most works in the literature are concerned with the stationary case and the recovery of 
isotropic sources / = f{x). Among those that address uniqueness and stability of reconstruc- 
tions in scattering media, we highlight [231 1211 [25]. We consider our work the counterpart 
of [23] for the time-dependent case. In fact, based on a Neumann series argument, analogue 
to that of [23], we have derived an iterative method for the recovery of the unknown initial 
condition. The method is convergent provided that the scattering coefficient is sufficiently 
small. This claim is made precise in Section |2] and then we prove it in Section [51 

For the stationary case, the use of the Neumann series is limited to certain smallness 
condition for the anisotropic portion of the scattering kernel as shown by Bal and Tamasan 
[25] . Stefanov and Uhlmann [21] by-pass this condition by reducing the inverse source 
problem to Fredholm form, showing that generic media an unknown isotropic source can be 
uniquely recovered in a stable manner from boundary measurements. Their proof assumes 
full boundary data and it is valid for all absorption and scattering coefficients in an open and 
dense subset of certain normed spaces. The case of partial data has been recently addressed 
by Hubenthal [25] whose approach is an extension of [24j. He is able to show that a source 
is recoverable if it is supported on certain regions that are visible from the accessible portion 
of the boundary. 

We also wish to mention some early works for special geometries or symmetries by Larsen 
[2n] and Siewert [27]. Larsen obtained a uniqueness result for sources embedded within a 
plane-parallel half-space. Siewert considered the transport problem in plane-parallel slab 
geometry and he obtained an algorithm to reconstruct the unknown source using spherical- 
harmonics expansions. A more recent practical work was carried out by Kim and Moscoso [9] 
for the problem in half-space with constant coefficients. Using appropriate Green's functions, 
they developed explicit formulae for the recovery of a point source or a piecewise constant 
source supported in a box. Most other studies for the stationary inverse source problem are 
concerned with non-scattering media and make use of the (attenuated) Radon transform and 
mathematical tools from integral geometry and microlocal analysis. An excellent literature 
review is found in [TUl Section 7]. 

1.2. The control problem. Control theory for PDEs is a very broad subject which has 
been investigated by a large number of researchers. We only direct the reader to following 
publications [2H1 [Ml [201 1311 [221 122] some which contain extended lists of references and 
overviews of important developments. From these references, it is clear that control theory 



has been comprehensibly developed for many PDEs of mathematical physics (wave, heat, 
Maxwell, elasticity, Schrodinger). However, this is not the case for the Boltzmann equation, 
not even in the linearized case. In fact, the exact controllability for transient radiative 
transport in heterogeneous media was first established in 2007 by Klibanov and Yamamoto 
|34] . We develop here an alternative proof for the exact controllability of the transport 
equation. As in the case of the inverse problem, we are particularly motivated by applications 
in medicine. For instance, radiation delivery for cancer therapy must be controlled in order 
to effectively destroy cancerous cell while minimizing damage to surrounding healthy tissues 
[35]. 

The exact controllability for radiative transport is to find inflow boundary conditions over 
a window of steering time to drive the transport field to a desired final state. The control 
data is usually not unique, but the first step is to establish its existence. The approach 
developed in [31] employs Carleman estimates which yield continuous observability even for 
time-dependent absorption and scattering coefficients, and also for optimal steering time. 
As a consequence, their results are much stronger than ours. However, we base our work on 
the analysis of the inverse problem, which is solved using time reversal. Therefore, our work 
reveals an alternative approach with added value worth reporting in this paper. The precise 
statement with regard to exact controllability is given in Section |2] and we provide a proof 
in Section [61 

2. Notation and statement of main results 

In this section we state the direct problem for transient radiative transport and its asso- 
ciated time-reversed problem. We also review some preliminary facts in order to state our 
main results in the proper mathematical framework. The formulation pursued here allows 
for heterogeneous media modeled by coefficients with low regularity. 

We assume that ^2 C M" is a bounded convex domain with smooth boundary dQ. The unit 
sphere in R" is denoted by S. The respective outflow and inflow portions of the boundary 
are 

{dn xs)± = {{x,e)ednxs : ± u{x) ■e>o} 

where u denotes the outward unit normal vector on dQ. It will be assumed that the particles 
travel at a flxed speed c > 0. In the rest of the paper, we will often make reference to the 
following scales: 

- / = diam(r2), the diameter of the bounded region Q, 

- c > 0, the speed of propagation, and 

- T = diam(r2)/c. 

Now we deflne the appropriate Hilbert spaces over which the radiative transport problem 
is well-posed. First, we denote by V'' and V^ the completion of C^{Q x S) with respect to 
the norms associated with the following inner products, 

(1) {u,w)yO = {u,w)L2^Qy,g) 

(2) (m,w)vi = 1^(0 ■ Vu,6 ■ Vw)yo + {u,w)yo + 1{W ■ 0\u,w)L^(^Qnxs) 

where V denotes the weak gradient with respect to the spatial variable x ^ Q. Now, denote 
by T the trace space deflned as the completion of C{dQ x S) with respect to the norm 
associated with the following inner product, 

(3) {u,w)T = l{\i^-9\u,w)L^(^dnxS)- 
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We also have the spaces T± denoting the restriction of functions in T to the in- and out-flow 
portions of the boundary dQ x §, respectively. Functions in V^ have well-defined traces on 
the space T as asserted by the following lemma whose proof is found in [3l IH |36l [371 ES [39] . 

Lemma 2.1. The trace mapping u i— )■ u\Qn defined for C^iVt x S) can he extended to a 
hounded operator 7 : V"*^ — ?■ T. Moreover, 7 : V^ — )■ T zs surjective. Analogous claims hold 
for the partial trace maps 7-1- : V^ —t- T±. 

In addition, we have the following definition for traceless closed subspaces of V^. 

(4) Yi = null{-f±) = {v eY^ : j±v = O} . 

For simplicity, we only pose and analyze the transport problem with vanishing incoming 
flow. In other words, we shall work in the space Vi. When a prescribed incoming flow 
needs to be included, it is easy to lift it as right-hand side source using the surjectivity of 
the trace operator 7_. The transient radiative transport problem for general heterogeneous, 
scattering media is the following. 

Definition 2.2 (Direct Problem). Given initial condition uq G Vl. and forcing term f G 
Ci([0,oo);V°) U C([0,oo);V^), find a solution u E Ci([0, 00); V^) n C([0,oo);Vi) to the 
following initial houndary value prohlem 

1 dv 1 

(5) -— + (9-V)u + fi^u + fiJl-K:)u=-f zn (0,00) X (fix §), 
c at I 

(6) u = uo on {t = 0} X (fi X §). 

We have included the factor l/l on the right-hand side of ([5]) so that both u and / have the 
same physical units. Here again V denotes the gradient with respect to the spatial variable 
X E fl. Here fi^ and /Xg are the absorption and scattering coefficients, respectively. The 
scattering operator /C : V° — )• V° is given by 

(7) {}Cu)ix,9)= f K{x,9,9')uix,e')dSi9'), 

Js 

where k is known as the scattering kernel. Throughout the paper we will make the following 
assumptions concerning the regularity of the absorption and scattering coefficients, and the 
scattering kernel. 

Assumption 2.3. The scattering coefficient < fig E L°°{Q). So there exists a positive 
constant /Zg such that < ^s{x) < /Is /o'^ ^-^^ ^ ^ ^- Similarly, the ahsorption coefficient 
< /ia G L°°{Q), with a constant Jl^ such that < /ia(2;) < /^a /^'^ ^-^^ x eVL. 

We consider a bounded and conservative scattering operator obtained by making the 
following assumption concerning the scattering kernel. 

Assumption 2.4. The scattering kernel < k G L'^{VL x S x §). It is also assumed that the 
scattering operator is conservative in the following sense. 



f k(x, 9, 9')dS{9') = 1, for a. a. {x, 9) E n x K 



In addition, we assume a reciprocity condition on the scattering kernel given hy 

(9) k{x, 9, 9') = k{x, -9', -9), for a.a. (x, 9,9') Eilx^xS. 

This means that the scattering events are reversihle in a local sense at each point x E fl. 



Now we turn our attention to a problem referred to as reversed transport. This problem will 
be used later in the analysis of the inverse problem. It is formally obtained by reversing both 
time t and direction 6 in the original transport equation ([5]), and by employing the reciprocity 
relation (^ to obtain the adjoint scattering operator /C*. This is expressed mathematically 
in dHJ. The reversed transport problem is defined as follows. 

Definition 2.5 (Time Reversed Problem). Given initial condition ipo G Y^ and forcing term 
p G C\[0, oo); V°) U C([0, oo); V^), find a solution ^ G C^i[0, oo); V°) n C([0, oo); V^) to the 
following initial boundary value problem 

Idih 1 

(10) -^ + (^ . V)^ - /iaV' - /i.(/ - /C*)^ = -p m(0,oo)x(fix§), 
c at I 

(11) ^ = ipo on {t = 0} X (fix §), 
where the adjoint scattering operator /C* : V° — t- V° is given by 

(12) (/C»(x,^)= f K{x,9',e)ijix,9')dSie'). 

Js 

Before presenting the main results of this paper, we wish to stress the importance of the 
reciprocity relation ([9]). This relation is derived from physical principles of scattering theory 
usually leading to rotationally invariant kernels of the form k = k{x, 9 ■ 9') which satisfy 
reciprocity. This relation plays a subtle but important role in this paper. It implies that 
the time reversed equation fITOl) (but not boundary conditions) coincides with the so-called 
adjoint equation. As a consequence, we obtain a simple but powerful relationship, expressed 
mathematically in fl56l) . which leads to the equivalence between the solvability of the inverse 
and the control problems. See the proof of theorem 12.91 found in Section [61 

2.1. Main result for the inverse problem. Now we state the inverse problem for transient 
transport along with its unique solvability and stabihty under the assumption that the 
scattering coefficient is relatively small. Our proof, presented in Section is based on a 
time reversal method inspired by the work of Stefanov and Uhlmann [40j. Our main goal is 
to provide a constructive proof that the initial state of the transport field can be uniquely 
reconstructed from time-resolved boundary measurements. For this inverse source problem, 
we assumed that the properties of the medium are known. 

Let u G C^([0,r]; V°) fl C([0,r]; V^) solve the direct transport problem 12.21 for unknown 
initial condition uq G N\. The outflowing boundary measurements are modeled by the 
operator A : Vi ^ C([0, r]; T+) defined as 

(13) (AMo)(t) = 7+wW, tG[0,r], 

where 7+ : V^ — )■ T+ is the out-flowing trace operator defined in lemma 12.11 With this 
notation we define the inverse problem as follows. 

Definition 2.6 (Inverse Problem). Let u be the solution to the direct problem \2.2\ for some 
unknown initial condition Uq and forcing term / = 0. The inverse source problem is, given 
the out-flowing measurement Auq, find the initial state Uq. 

Our main result concerning this inverse problem is the following. 

Theorem 2.7. Assume that //Ig e'^^=''^'^=-' < e~^. Then there exists a time r < 00 such that 
the out-flowing boundary measurement Kuq G C([0,r];T+) determines the initial condition 
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Uq G Y^ uniquely. Moreover, the following stability estimate holds, 

ll^oilvi < C'||A-Uo||c([0,r];T+) 

for some positive constant C = CijI^,]!^, I). 

This theorem is a consequence of theorem 15.11 which is stated and proved in Section |5l 
See also the remarks following the proof of theorem 15.11 concerning the validity of theorem 
12.71 in a weaker setting, that is, for initial data uq G V°. Notice that theorem 15.11 provides 
a convergent iterative method for the reconstruction of the unknown initial condition uq. It 
is also worth mentioning that the stability estimate of theorem 12.71 is optimal for the chosen 
norms because we can easily show that ||mo||vi dominates ||Ano||c([o,r];T+) up to a constant. 
In other words, in this stability relation we cannot relax the norm of the measured data Auq 
without relaxing the norm of the initial condition uq as well. 

2.2. Main result for the control problem. Here we proceed to define the control problem 
and state our main result concerning exact controllability of the transport field from control 
boundary data. Our proof, based on duality arguments, is presented in Section O 

For the controllability issue, we work within the framework of a Hilbert space for the 
control functions on the boundary (generalized traces) and the corresponding mild solutions 
of the radiative transport problem. For the existence of mild solutions in semigroup theory, 
see the standard references [IT] H2] . The treatment of generalized traces for mild solutions 
can be found in |3ll Section 2] or [H Section 14.4] which is based on Cessenat [371 EH]. See 
also Bardos US pp. 205-208]. 

We consider the following transport problem with prescribed infiow data. Given h G 
L^([0, r]; T_), find a mild solution v G C([0, r]; V°) of the following problem 

1 dv 

(14) -— + {9-V)v + fi,v + fi,{I-}C)v = in(0,r] X (fixS), 

(15) v = on {t = 0} X (fi X §), 

(16) v = h on [0,r] X (9fi X §)_. 

The infiowing boundary control is modeled by the bounded operator T : L^([0, r]; T_) — )■ V° 
defined as 

(17) Th = v{T), 

where v is the mild solution of the problem flT^ - flT6l) . With this notation, we define the 
control problem in precise terms as follows. 

Definition 2.8 (Exact Controllability). Given a target state f^ G V°, find a finite steering 
time r > and an inflow control condition h G -^^([0, r]; T_) such that Th = v^. 

Our main result concerning this control problem is the following. 

Theorem 2.9. Assume that Ijl^e^^^'^^'^^^ < e^^ . Then there exists a steering time r < oo 
such that for a given target state f^ G V°, there exists inflow control h G L^([0, r]; T_) so 
that the mild solution v G C([0,r]; V°) of the problem ([77])-( l7^) satisfies v{t) = v^. 

The control boundary condition has the form h = /imin + 9 where g G null(T) and /imin G 
null(T) is uniquely determined by f* as the minimum-norm control satisfying, 

||^mm||L2([0,T];T_) < C'||'?^*||vO 

for some positive constant C = C(jI^,JIg, I). 



3. Analysis of the direct problem 

In this section we briefly review the well-posedness of the direct problem 12.21 and the 
reversed problem 12.51 in appropriate functional spaces. Most of the results in this section are 
well-known or easily derived from the literature for radiative transport equations. However, 
we explicitly state the key ideas in order to use them later in the analysis of the inverse and 
control problems. Our approach is mainly based on semigroup theory and its application 
to evolution PDEs [H HH |12]. Besides reviewing some standard results, the purpose of 
this section is to specify conditions on the absorbing and scattering coefficients that make 
the semigroups associated with the two evolution problems decay exponentially fast. See 
proposition 13.51 at the end of this section. It will become clear in the analysis of the inverse 
problem (see Section [5]) that the exponential decay of these semigroups is essential for our 
proof of the main results of this paper. 

In order to accomplish our goal for this section, we will treat the transport problems as 
bounded perturbations of the ballistic portion of the transport equation. Hence, we start by 
analyzing the ballistic (non-scattering) transport problems. 

Definition 3.1 (Ballistic Problems). For the direct problem, given initial condition Uq G Vi, 
find a solution u G C^([0, oo); V°) fl C([0, oo); VlL) to the following initial boundary value 
problem 

1 du 

(18) -— + (0- V)M+(/ia + /is)M = m(0,CX)) X (fixS), 

c at 

(19) u = uq on {t = 0} X (fi X §). 

Similarly, for the time-reversed ballistic problem, given initial condition ipo € V^, find a 
solution ip G C"'^([0, cxd); V°) fl C([0, cxd); V^) to the following initial boundary value problem 

(20) -^ + (^-V)^-(/ia + /is)V' = zn(0,oo)x(fix§), 
c at 

(21) ^ = ^0 on {t = 0} X (fi X §). 

The ballistic problems are well-posed. In fact, solutions can be written explicitly using 
the method of characteristics. Since we plan to use the theory of semigroups, then we pose 
these ballistic problems as abstract Cauchy problems in Vi with the following operators 
Aq,Bq:Y\^ V° given by 

(22) A^v = -c [{e ■ V)v + (/ia + f^s)v] , 

(23) Bov = -c [{e ■ V)v - (/ia + fis)v] , 

as their respective generators. We have the following properties concerning the strongly 
continuous semigroups generated by Aq and Bq. 

Lemma 3.2. Let {So{t)}t>o and {-Ro(^)}t>o be the Co-semigroups generated by Aq and Bo, 
respectively. Then the following properties are satisfied. 

(i) Let V G Y^ be fixed, then the maps t H- So(t)v and t i— )■ Ro{t)v are continuous from 

W^ to Vi and continuously differentiable from W^ to V°. 
(ii) For each t > 0, both So{t) and Ro{t) extend as bounded operators from V° to V°. 
(iii) For all t>0, we have \\So{t)\\ < 1 and \\Ro{t)\\ < e'(^-+'^=). 
(iv) For all t > T we have So(t) = Ro{t) = 0. 
(v) For allt>0 and allu < 0, we have \\So{t)\\ < e'^^*"^) and \\Ro{t)\\ < e^(*-^)e'(Pa+7is). 
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Proof. The proof of properties (i)-(iii) follow from the standard theory of semigroups [^H2] 
and the fact that the solutions of the Cauchy problems can be explicitly expressed using 
the method of characteristics. Similarly, property (iv) holds because we impose vanishing 
inflow condition on the boundary of Q by working in the space Vi. Thus, for time t > T 
the support of solution u = u{t) or ip = ip{t) has already left the domain Q. Finally, the 
combination of (iii)-(iv) yields the estimates in (v) as desired. D 

Now we turn to the full problem in heterogeneous, absorbing, scattering media. We view 
both, the direct problem 12.21 and the reversed problem 12.51 as perturbations of the ballistic 
problems 13.11 So we pose the former two problems as abstract Cauchy problems in V^. 

The direct problem 12.21 corresponds to 

(24) u{t) = Au{t) + ^ /(t) for t > 0, 

(25) «(0) = Mo, 

where T = l/c and the operator A : Vi C V*^ — )■ V° is given by 

(26) Av = -c[{e-V)v + fi^v + i2s{I-IC)v]. 
Similarly, the reversed problem 12.51 corresponds to 

(27) 7p{t) = B^Pit) + ^ p{t) fort>0, 

(28) V(0) = V-o, 
where the operator B : Vi C V*^ — )■ V° is given by 

(29) Bv = -c [{9 ■ V)v - fx^v - fi,{I - )C*)v] . 

In order to apply the theory of bounded perturbation of semigroups [IT, Ch. 3], we find 
it convenient to state and prove the following lemma. 

Lemma 3.3. Let /C, /C* : V° — )► V° he defined by ^ and IH^) . respectively. Then we have 

II iK"ll II V* II 1 

||A^|| = ||A^ II = i. 

Proof. Let n G V° and consider, 

I|/Cm||^o = J^ /s (/s k{x, 0, e')u{x, e')de'fdedx 

< j^ /g /g k{x, e, d')dd' /g k{x, 0, 0')u{x, e'fde' dodx 
= J^ /g /g k{x, e, e')u{x, e'fde'dOdx 

= J^ J^ Jg k{x, 9 , 9')d9 u{x, 6'yd6'dx = \\u\\lo- 

Here we have used the Cauchy-Schwarz inequality, Fubini's theorem and assumption 12.41 on 
the scattering kernel k. Now, notice that equality is attained for the function u{x,9) = 1. 
Therefore, it follows that ||/C|| = ||/C*|| = 1. D 

Now we may state and prove the well-posedness of the direct problem 12.21 and reversed 
problem 12.51 in terms of the associated Cauchy problems fl25|) - fl26l) and (!28|) - (!29|) . respectively. 

Theorem 3.4. The direct problem\EM has a unique solution given by the formula 

1 /■* 

u{t) = S{t)uo + - S{t-s)f{s)ds, 



where {S{t)}t>o is the Co-semigroup generated by A satisfying properties (i)-(ii) of lemma 
[23 and such that S{t) : V° ^ V° zs bounded with \\S{t)\\ < M^e('^+*^""^=)* for all t > and 
allcj <0 where M^ = e"'^^. 

Similarly, the reversed problem \2.5\ has a unique solution given by the formula 

ij{t) = R{t)iJo + - / R{t- s)p{s)ds, 

where {R{t)}t>o is the Co-semigroup generated by B satisfying properties (i)-(ii) of lemma 
[23 and such that R{t) : V° ^ V° is bounded with \\R{t)\\ < AT^e^'^+^-'P^)* for all t > and 
alloj <0 where N^ = e'^^T ^K-p^+ji^) ^ 

Proof. The proof for both problems 12.21 and 12.51 is exactly the same, so we only address the 
well-posedness of the direct problem 12.21 We use [HI Ch. 3 : Thm. 1.3] and property 
(v) in lemma 13.21 to obtain the existence of S{t) as a Co-semigroup satisfying properties 
(i)-(ii) of lemma 13.21 and the estimate on its norm ||S'(t)||. Now since we assume that 
/ e C([0,oo);Vi) UCi([0,oo);V°), the regularity of the term 

I S{t-s)f{s)ds, 
Jo 

follows from [41, Ch. 7: Section 7]. D 

The estimate ||-R(t)|| < A^^e*^'^+^'^'^'^=^* found in theorem 13.41 motivates us to find conditions 
on 71^, /Ig and a good choice of a; < such that the reversed semigroup R{t) becomes a 
contraction for sufficiently large time t > 0. It will be evident that such a property determines 
the solvability of the inverse problem. See below in Section [5] So consider the function 
i? : R — )■ M given by the exponential rate found in the estimate ||-R(t)|| < N^e^'^^^"'^'^^'^^ , or 
in other words, 



E(u) = u + e-'^^e'('^-+^=) 



CfX, 



s" 



We would like to find an optimal choice u* < that minimizes E over R. Subsequently, we 
would like to find a condition on fl^ and fX^ leading to E(uj*) < 0. Clearly, u* is the unique 
solution to E'{u) = which is easily obtained to be 



u 



In (Te'(^-+^=)c/i,) 



T 
Recalling that / = To, it follows that 

E(uj*) = — ^ ^-^ + - = — ^ ^-^-^ 

\ / rri ' rri rri 

Hence, in order for E{u*) < 0, we require that Ijl^e^'^^'^'^^^'^ < e^^. Now, from the estimates 
in theorem 13.41 we clearly arrive at the following result. 

Proposition 3.5. Let I = diam(i7). // 

/7x,e'(^-+^=) <e-^ 

then both ||>S'(t)|| and ||-R(t)|| decay exponentially fast according to 

\\S{t)\\ < \\R{t)\\ < ee'(^-+'^=) (ee'('^-+^=) IJ^sY^'''^ , t>T. 
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4. Analysis of the stationary problem 

In this section we state the stationary or steady-state problem for the transport equation in 
general heterogeneous, absorbing, scattering media. This is done for the stationary problems 
corresponding to both, the direct problem 12.21 and also the reversed problem 12.51 We shall 
prove the well-posedness of both problems under the same assumptions of proposition 13.51 

In the analysis of the inverse problem, it will become clear that the well-posedness of the 
reversed stationary problem plays two important roles. The first of these roles has to do with 
the definition of a time reversal operator that solves the inverse problem up to a contraction 
map. More precisely, this operator is written in terms of a time-reversed evolution problem 
which in turn needs an initial state to be well-defined. The needed initial state is then 
provided by the solution of a reversed stationary problem in order to obtain a crucial stability 
property. The second role is that a carefully chosen solution for a reversed stationary problem 
will allow us to stay within the framework of the space V^. As a consequence, the stability 
of the source reconstruction method is given in terms of the V^-norm which is stronger than 
that of the space V°. 

The reversed problem is particularly challenging since the coefficients in the PDE appear 
with the reversed signs. This means that the term containing /la acts as emission instead of 
absorption. As a result, the weak formulation of this problem does not lead to a coercive 
(uniformly convex) form. Instead, it leads to an indefinite (saddle-point) problem. We 
shall overcome this difficulty by employing the celebrated theorem of Babuska [H] . This is a 
generalization of the Lax-Milgram lemma especially designed for such saddle-point problems. 
Our approach follows the ideas developed in [5SlEn], but we introduce certain simplifications 
and modifications. 

Now we proceed to state the direct stationary problem with vanishing incoming fiow. In 
other words, we shall work in the space Vi. As in the case of transient transport, when a 
prescribed incoming fiow needs to be included, it can be lifted using the surjectivity of the 
trace operator 7_. 

Definition 4.1 (Direct Stationary Problem). Given / G V° find a solution u G Y\ to the 
following boundary value problem 

(30) {e-V)u + fi^u + fi,{I-}C)u = -f mfixS. 

And the reversed stationary problem is defined as follows. 

Definition 4.2 (Reversed Stationary Problem). Given p E Y^ find a solution ip G V^ to 
the following boundary value problem 

(31) [9 ■ V)^ - /iaV' - yUs(/ - /C*)V' = yp m fi X §. 

One of the most important tools in our analysis is the following integration-by-parts for- 
mula or Green's identity. For functions m, f G V^ we have 



(32) / {e-Vu)v= / {d-u)uv~ / {d-Vv)u, 

JqxS JdQxS JnxS 

and in particular if we let f = n G Vi and multiply by /, we obtain from Young's inequality 
that 

(33) ^11(0 . V)u\\lo + h\u\\lo > f lie- Vu)u = U (0. uy = hu\\l 
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We will also need a stability estimate for the ballistic portion of the stationary problems. 
This may be obtained from a Poincare-type inequality to estimate the V°-norm of a function 
in terms of its directional derivative provided that inflow values vanish. Since the ballis- 
tic solution can be written explicitly using an integrating factor, the proof of a Poincare 
inequality from |39] can be easily modified to obtain the following stability estimates. 

Lemma 4.3. Let I = diam(r2). For all v G Y^ we have that 

||f llvo < \/2l \\{6 ■ V)f + (/ia + /^s)f llvo 

||(^-V)^;||vo < (^l + y2/(7I, + 7lj) \\ie-V)v + ifi^ + fis)v\\vo 

and 

\\v\\yo < 72/ e^(''-+^=) \\{e-V)v-{i2, + ij,)v\\Y0 
||(^-V).;||vo< (l + v^/(/I, + 7lJe'(^^+^=)) ||(e.V)t;-(/ia + /is)t;||vo- 

In order to prove the well-posedness of the stationary problems 14.11 and 14. 2[ we first set 
up associated variational problems similar to those of [38]. However, as opposed to [38j, we 
directly seek for a solution in the trial space Vi which already has enough regularity for the 
solution to satisfy the transport PDE in a strong sense. Without further ado, we define the 
bilinear forms governing the stationary problems I4.HI4.2[ 

Let a, b : Y^ x V° — )■ M be a bilinear forms given by 

(34) a{u,v) = I {{9 ■ 'V)u + fiaU + fis{I — JC)u,v)yo, 

(35) 6(^,0) =/((e-V)^-/ia^-/is(/-/C*)^,0)vo. 

Hence, the stationary problems 14.1114.21 are equivalent to find u^ip E N\ such that 

(36) a{u,v) = {f,v)YO for all v e V°, 

(37) bitfj, (f)) = (p, 0)vo for all G V^. 



As stated above, the reversed problem 14. 2[ governed by the bilinear form 6(-, ■) poses the 
greater challenge due to its lack of positive definiteness. Hence, we first proceed to analyze 
this case, and the steps can be easily modified to deal with the more favorable structure of 
the other bilinear form a(-, ■). We start with the following lemma. 

Lemma 4.4. // //I^ e'*^^"^^^-* < e~^ , then for each ip G Y^ there exists </> G V such that 
b{4',(p) > /5||^||vi||0l|vO; where the constant (3 > is independent of ip and (p. In fact, 



2v^ max{/3i,/32}^ 



where /3i = (l + V21{JI.^ + /Ij e'(^-+^=)) and P2 = V^ €^^''-^''-1 

Moreover, for each non-zero G V° there exists ip G Vi such that h{ip^ 0) > 0. 

Proof. For given ip G Vi, let = / {{6 ■ V)ip — (/ia + Ps)^) ^ V°. Then we have, 

6(7/^,0) = 11011^0 +/(Ps/C*V',0)vO 

>ll*-//lj|^||vo||0||vo 

> (1 - y2/7I,e'(P^+P=)) 11011^0 > (1 - v^/e) 
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where we have used an estimate from lemma I4.3[ Now it only remains to show that 
dominates HV^Hvi- From the definitions of (3i and P2 and lemma H73| we obtain that 

||0||vo > U\iO ■ V)^||vo and ||0||vo > -^l|^||vo. 

Pi P2 

Finally, using Green's identity ( 133|) we arrive at the following inequality 

The desired result follows by recalling the definition of the norm in V^ given in ([2]). 

For the second part, given non-zero G V°, we choose ip € Yl_ such that {6 ■ V)ip — (/ia + 
/is) = 4>/l in (fi X S). This can indeed be accomplished since the ballistic stationary problem 
with vanishing infiow data is uniquely solvable in Vi. Therefore, we have 

K^A) > 11011^0 -/7lj|^||vo||0||vo 

> 11011^0 (1-v^ //Is e'(^^+^=)) 
>||0||2o(l-v^/e)>O, 

where we have used an estimate from lemma 14.31 This concludes the proof. D 

The above lemma holds under the condition that //Ige'*^'^""*"^^-* < e~^. This choice was 
purposely made to coincide with the hypothesis of proposition 13.51 It will become clear in 
Section [5] that the conclusion of proposition [375] is crucial for the well-posedness of the inverse 
problem. 

Lemma [4.41 and Babuska's theorem [33] lead to the well-posedness of the reversed station- 
ary problem 14.2! as stated in the following theorem. 

Theorem 4.5. If 1]!^ e^-'^^'''^^-' < e~^, then for each p G V°, the reversed stationary problem 
4.^ has a unique solution ip G Vi such that / f57]) is satisfied a.e. in (i7 x S). Moreover, the 



solution satisfies the following stability estimate, 

IIV'llvi < ^IIpIIvo, 



where (3 > is defined in lemma\4^\ Same conclusion holds for the direct stationary problem 



4.1. 



This theorem implies that the so-called subcriticality condition, /ia(a;) > 0, assumed in the 
literature is not necessary for the stationary radiative transport problem to be well-posed. 
Notice that the proof of lemma 14.41 is indifferent to the sign of the coefficients /ia and /is, 
provided that Ijig e'*-^""'"^''-* is sufficiently small. This last remark implies the well-posedness of 
the direct stationary problem |4. II for Ifl^ ^^(Ma+'Us) < g^^. However, other proofs can be found 
in the literature since the direct stationary problem 14.11 satisfies the subcriticality condition. 

We also wish to state here the following corollary that we will employ in the analysis of 
the inverse problem in Section [5] 

Corollary 4.6. Let {S{t)}t>o o-nd {R{t)}t>o be the Co-semigroups associated with the direct 
and reversed problems \2.2\ and \2.5[ respectively. Ifljl^e^^"-^'^'''' < e^^, then these semigroups 
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are bounded linear operators from V^ to Y^ satisfying the following estimates 

\\S{t)\W. < i+^yP.) ||g(t)||^o and \\R{t)\W^ <^-^^^^^f^^^\\R{t)\Wo 



where P > is defined in lemma \4^ . 



Proof. Let Uq G Y^ be arbitrary and consider, 

ll'5(t)Mo|ki < I^||M(t)||vo = ^||AM(t)||vo 

= ^\\AS{t)uoho = j^\\S{t)Auoho 

< ^||^(t)||vopMo||vo < ^\\S{t)Uo\\A\\\\uo\\yi. 

We have also used the fact that a semigroup commutes with its generator when acting on the 
domain of the generator. In other words, S{t)Av = AS{t)v for all v G Y^. The boundedness 
of S{t) : V*^ — )■ V° was obtained in theorem 13.41 and it is clear that A : Vi — ^ V° is bounded 
with ||y4|| < c(/~^ + 7^a + 27is). The proof for R{t) is similar. D 

5. The inverse problem 

With the tools developed in the previous two sections, we are ready to solve the inverse 
problem 12.61 We base our analysis on the time reversal method, which is usually employed 
for the wave equation. In particular, our approach mimics that of |10] for the thermoa- 
coustic tomography problem. The time reversal method is particularly well-suited for the 
wave equation because that equation is invariant under time inversion. Unfortunately, the 
radiative transport equation does not enjoy such a property; and this represents our main 
challenge. However, the theory expanded in Sections |3]|4] provides the necessary tools to deal 
with this difficulty. 

Let u G C\[0, r]; V°) n C{[0, t];YI) solve the Cauchy problem, 

(38) u{t) = Au{t) for < t < r, 

(39) u{0) = uo, 

where the operator A : Vi — )■ V° is defined in fl26|) and the initial condition uq G Vi. The 
boundary measurements are modeled by the operator A : Yl_ — )■ C([0, r]; T+) defined by (TT3|) 
for some chosen time r > 0. Recall also Definition 12.61 of the inverse problem. 

We wish to determine the existence of a finite measurement time r which ensures the 
recovery of Uq from Auq. Clearly we need r >T = diam(n)/c because T is the time required 
for the special case when the medium is non-scattering. 

In order to employ a time reversal argument, we define the following reflection operator 
U : C([0,r]; V^) -^ C([0,r]; V^) given by 

(40) {Uv){x,e,t) = v{x,-e,T-t). 

This operator is certainly bounded, and in fact it is norm-preserving. Moreover, it is also well- 
defined and norm-preserving as U : Ci([0,r];V°) -^ Ci([0,r];V°) and U : C([0,r];T±) -^ 
C([0, r]; Tip), and in all cases U"^ = I. We will also use the fact that 

(41) UIC = ICU 

which follows from the reciprocity condition ([H]) satisfied by the scattering kernel k. 
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It is not hard to see that u := (Uu) G C^([0, r]; V°) H C([0,r]; V^) solves the following 
reversed initial boundary value problem, 

(42) u{t) = Bu{t) for < t < r, 

(43) ^(0) = (f/n)(0), 

(44) 7„M(t) = (f/AMo)(t) forO<t<r. 

where the generator B : Y^ — )■ Y^ is defined in (1291) . 

If we had access to {Uu){0), then we could solve (H21)- (H^ and it would follow that 
Mo = {Uu){0). Unfortunately, this is not realistic and we only have access to the boundary 
measurements Auq. For the time reversal method, {Uu){0) is simply replaced by a known 
function tp of our choice. 

Inspired by the work of Stefanov and Uhlmann [ID] we employ the reversed stationary 
problem 14.21 as a lift to obtain an initial condition ip that conforms to the boundary data 
{UAuo){t) at t = 0. Such a function would satisfy a boundary value problem of the form, 

(45) (^ ■ V)^ - /iaV' - yUs(/ - /C*)^ = in(nx§), 

(46) j_ip = ho on (dn x §)_. 

for ho eT_. As seen below, in practice we choose ho = {UAuo){0). 

With this choice of ip, we proceed to define the time reversal operator which acts as an 
approximate left-inverse for A. Given h e C([0, r]; T+), find v e C\[0, r]; V°) H C([0, r]; V^) 
satisfying 

(47) v{t) = Bv{t) for < t < r, 

(48) v{0) = ^, 

(49) ^_v{t) = {Uh){t) forO<t<r. 

where ip satisfies (05])- (06]) with ho = {Uh){0). 

The time reversal operator G : C([0, r]; ¥+) — )■ V"*^ is given by 

(50) Gh={Uv){0). 

We should prove the well-posedness of the boundary value problem (l45!) - (H6l) and the 
initial boundary value problem (H7|) -f H9|) . These problems are not the same as the respective 
reversed problems 14.21 and 12.51 because now we are prescribing a non-zero inflow data in P6ll 
and (149|) . However, since the trace operator 7_ is surjective, we can always lift the boundary 
data and pose new problems in the Y^. Then, theorems 14.51 and 13.41 yield the well-posedness 
of these two problems, respectively. As a consequence, then G is a bounded operator. 

We are interested in making I — GA a contraction mapping in the space Vi for some 
properly chosen measurement time r < oo. This is the basic idea employed in |1D]. We 
proceed to state our result in the form of a theorem. 

Theorem 5.1. Assume that Ifl^ e^'Pa+'Us) < g-i, i^t Q ; y^ —^ Y^ be given by Q = I — GA. 
There exists a final time r < oo such that Q : Y^ — t- Y^ is a contraction and {I — Q) : Y\_ — > 
V^ is boundedly invertible. Moreover, the solution to the inverse problem \2.6[ is given by 



Uo 



Y^ Q'^Gh, h = Auo, 



n=0 



with convergence in the Y^-norm. 
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Proof. Following the approach in gOl Thm. 1], let w G C\[0,t];¥^) n C([0,r]; Vi) solve 
the following Cauchy problem, 

(51) w{t) = Bw{t) for < t < r, 

(52) w{0) = (f/M)(0)-^G Vi. 

Notice that this is in the precise form of the reversed Cauchy problem 12.51 and the initial 
condition in ( 152|) truly belongs to Y\ by design. Hence, theorem 13.41 implies that w{t) = 
R(t){{Uu){0) — ip). Since ( 15T1) is satisfied in a strong sense at t = r, then it follows from 
theorem 14.51 that 

(53) ||^(r)||^. <l|mr)||vo. 

cp 

Notice also that {v + w) and {Uu) solve the same initial boundary value problem (I42l)-( l44l) . 
By uniqueness then w = Uu — v for all t G [0, r], in particular if we apply U and evaluate 
at t = we obtain (f/w)(0) = Uq — GAuq. We have used the definition of G given in (!50|) . 
Then we obtain that Quq = {Uw){0). Therefore, the following estimates hold 

lignollvi = \\{Uw){0)\Wi = \\w{T)\\yi < ^||w;(r)||vo = ^\\Bw{T)\\yO 

= ^||5i?(r)((f/n)(0) -^)||vo = ^||i?(r)i?((f/n)(0) -^)||vo 

= i-^\\R{T)B{Uum\Wo < ^||i?(r)||vo||i?||||(f/w)(0)||vi 

= ^^\\R{r)\Wo\\B\\\\u{r)\W^ = ^||i?(r)||vo||i?|| ||5(r)«o||vi 

<^p(r)||vo||i?||^||S(r)||voP||||wo||vi. 

We have used the norm-preserving properties of the reflection U. The estimate 053p was used 
in the second inequality and the constant /3 is defined in lemma W^ Recall that Bip = from 
fH5|) . We have also used the fact that a semigroup commutes with its generator when acting 
on the domain of the generator. In other words, R{t)Bv = BR{t)v for all v G Y^. The 
boundedness of S{t) and R{t) was obtained in theorem 13.41 as maps on V°, and in Corollary 
14.61 as maps on V^ . Finally, it is clear that A :Y^ ^-Y^ and B :Y^ ^ Y^ are bounded with 
II ^11 = ||-B|| < c{l~^ +JI^ + 2/1^). Hence we obtain the following estimate, 

WQuoh^ < (llIit^^tltA^ ||i?(r)||vo||^(r)||vo||«o||vi 

Now the stability estimate from proposition 13.51 implies that both ||S'(r)||vo and ||i?(r)||vo 
decay exponentially fast as r increases. Therefore, there exists a finite time r such that 
Q : Yl_ — )■ Y]_ is a contraction mapping and the desired results follow from the Neumann 
series theorem. D 

We point that theorem 15.11 directly implies the validity of the first main result expressed 
as theorem 12.71 Notice that the choice of ip in ( HHI) as the solution of the reverse stationary 
problem fHSl) - fH^ was crucial in order to stay within the formulation of the space V^ and 
to obtain the estimate on the norm of Q. These two facts follow from the following two 
properties: 

- i/j conforms to the boundary data (f/Auo)(t) at t = 0, and 

- ip belongs to the null space of the generator B. 
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From the well-posedness of the reversed stationary problem \A.2\ we see that the above two 
properties determine ip uniquely within the space V^. 

Now to conclude this section, we wish to mention that the time reversal operator G 
from ( 15U]1 can be modified to obtain a bounded operator G : L^([0, r]; T_|.) — )■ V°. This is 
accomplished using the concepts of mild solutions and generalized traces as in [ST, Section 
2], [U Section 14.4] or Cessenat [3Z1 ES]- In fact, the measurement operator A from (TT51) 
can be boundedly extended to A : V° — )■ L^([0, r]; T_|_) as shown in the proof of theorem 12.91 
in Section [61 However, it would no longer make sense to speak of measured boundary data 
{UAuQ){t) at time t = 0. Hence, the function ip and the operator G have to be modified 
as follows. We let v e C([0,r];V°) be the mild solution of (I47])-(149]) with ^ = and 
h G L^([0, r]; T+). Then, we can obtain (following steps analogous to those in the proof of 
theorem 15. ip an estimate on the norm of Q : V*^ — ?■ V°, 

IIQuollvo < ||-R(T")||vo||'S'(r)||vo||uo||vo, for all Uq e V°. 

In view of proposition 13. 5[ we see that there exists r < oo such that Q : V^ — ?■ V is 
a contraction mapping and the conclusion of theorem 15.11 is valid with convergence in the 
yO-norm for initial condition uq G V°. Finally, we also have that the conclusion of theorem 
12.71 is also valid in a generalized sense, meaning that 

(54) \\uo\\yo < C||AMo||L2([o,r];T+) for all Uq e V°, 
for some positive constant G = G(jI^,]Ig,l). 

6. The control problem 

In this section, we develop the proof of our main result concerning the exact controllability 
problem the radiative transport. It is well-known in control theory that exact controllability 
in a Hilbert space setting is equivalent to the continuous observability property for the so- 
called adjoint problem. In turn, observability is obtained from the solvability of the inverse 
problem which we have already established in Section [51 The relation between these concepts 
is made precise using duality arguments which we proceed to describe. 

We will need the following angular-reflection operator \^ : V° — t- V° defined by 

(55) {Vw){x,e)=w{x,-e), 
which is clearly unitary. 

Proof of theorem \2.9i . We wish to construct the adjoint of the measurement operator A de- 
fined in (iT3l) viewed for now as a densely defined operator A : V^ C V° — )■ L^([0, r]; T+). 
First, let Uq G Vi be arbitrary and u G G^{[0,t];Y^) fl C([0,r]; Vi) be the unique solution 
of the problem ([38l) - (!39l) with uq as initial condition. Hence we have that 7+M = Auq. 

Now let h G L^([0,r];T+) be arbitrary, and {/ifc}fc>i C C^([0, r]; T+) be a sequence con- 
verging to h in the L^{[0, r]; T+)-norm. For each hk let Vk G G\[0, rj; V°)nC([0, r]; V^) be the 
strong solution of the problem (fT4l) -( !T6l) . Hence, Vkir) = Th^ where T : -^^([0, r]; T+) — ?■ V° 
is bounded as defined in ( ITTll . 

Recall the reflector operator from ( HUj) and notice that = Uu solves 



0(t) = 50(t) for < t < r, 

0(r) = {Uu){t), 

-f+(f){t) =0 for < t < r. 
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Now integrate = {vk — Avk)4> over the domain [0, r] x Q x S and use integration-by-parts 

to obtain (wfc(r), 0(r))vo - (t;fc(0), 0(O))vo = (t^fc, 0)l2([o,t];T_) - (^^fc, 0)l2([o,t];T+)- Recall that 
'^fcl'?") = '^hk, Vk{0) = 0, 7_ffc = /ifc, 7+0 = 0, and 7_0 = 'y^Uu = U'y+u = UAuq. Also 
notice that 0(r) = Vmq- It follows that {Thk,Vuo)YO = {hk, UAuq) L2(^[Q^^yj_y Recall that T 
is bounded, so after taking the limit A; — ;■ oo we obtain 

{Th, Vuo)yo = {h, f/Auo)L2([o,T];T_); ^oi all h e L'^{[0, r]; T+) and Uq E Vi 

where U and V are the reflector operators defined in ( HUj) and ( 15^ . respectively. Because 
V^ = V"* = 1/-1 and U = U* = U-\ we find that 

(56) A* = VTU 

Since V^, T and U are bounded operators then so is A* and consequently A can be boundedly 
extended to A : V° — ;■ L^([0, r]; T+). We have already used this extension in Section E] to 
obtain f l54p . This means that A : V° — )■ -^^([0, r]; T+) is injective and it has a closed range. 
Because T* = UAV, with U and V being boundedly invertible, it follows from the closed 
range theorem that both range{T*) C L^([0, r]; T_) and range{T) C V° are closed. 

The exact controllability problem 12.81 reduces to showing the surjectivity of the control 
operator T. Since range(T) is closed in V°, then basic duality theory tells us that range{T) = 
null{T*) = null{UAV) . Now, we know that U and V are isometrics and A is injective. 
Hence range{T) = V°, which establishes the exact boundary controllability of the transport 
field. Moreover, from pseudo- inverse theory for Hilbert spaces [15], we obtain the minimum- 
norm control given by h^^^ = T*(TT*)~^f^, which concludes the proof. D 
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